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Abstract 

We review the current status of D°-D° mixing, with 
special emphasis in the most recent results. We begin 
with a discussion of charm mixing and CP violation phe- 
nomenology, the evolution with the decay proper time, 
and physics processes contributing to these. Then we 
follow with the summary of the main experimental tech- 
niques and the results in the various final states. We 
make use of the analysis reporting the first evidence of 
D°-D° mixing by BaBar in D° — > K + tt~ decays as a 
textbook example, and then we discuss the results from 
other two-body and three-body final states. We conclude 
with the combination of all experimental results. Time- 
integrated CP violation measurements are not discussed 
here. 

1. Introduction 

Particlc-antiparticlc oscillation (also referred as mix- 
ing) is a well known phenomenon observed in the kaon 
system in 1956 [Tj, in the B% system in 1987 [2J, and more 
recently in 2006 in the B® system [3J. Mixing and CP 
violation (CPV) in the charm sector were first discussed 
over three decades ago [3], but experimental evidence 
for oscillation has been presented only in the last two 
years [SJ H El [HI IS HD] , and no evidence for CPV has yet 
been reported, with upper limits currently at about 1% 
level. 

Charm mixing is the only involving down-type quarks 
in the mixing loop, since neutral pions do not oscillate 
and the top quark does not have bound states. Thus 
D°-D° mixing offers an unique probe for New Physics 
(NP) via flavor changing neutral currents (FCNC) in the 
down-quark sector, providing interesting constraints on 
NP models. The caveat is how to distinguish NP from 
Standard Model (SM) long-distance (non-perturbative) 
uncertainties. A possible avenue is correlating charm 
mixing studies (and possibly also rare charm decays) 
with a comprehensive account of CP violation in D°- 
D° mixing (both within the SM and beyond). The two 
effects are heavily suppressed in the SM (charm mix- 
ing is about two orders of magnitude slower than in 
the neutral-kaon system and CPV is well below the 
per mille level), which makes these experimentally dif- 
ficult to observe, although NP can produce significant 
enhancements. 

2. Charm mixing phenomenology 

Neutral- D mesons are created as flavor eigenstates of 
strong interactions, but they mix through weak inter- 
actions. The time evolution is obtained by solving the 
time-dependent Schrodingcr equation, 



with H w = M— iT/2 the effective Hamiltonian, where M 
and r are 2x2 matrices that represent transitions via off- 
shell (dispersive) and on-shell (absorptive) intermediate 
states, respectively. Assuming CPT invariance, we have 
Mix = M22 and Tn = r 2 2- Since these matrices are 
Hermitian, M 12 = M 2 i and Ti 2 = T^. If CP in mixing 
is conserved, Mi 2 = M21 and Ti2 = T2i. 

The physical (mass) eigenstates are linear combina- 
tions of the interaction eigenstates, |X?i 2) = p\D°) ± 
q\D°) [IT], with time evolution \D h2 (t)} = e~ lXl - 2t \D h2 ), 
where Ai,2 = m.1,2 — ^1,2/2 are the eigenvalues. Here, 
m i,2 (ri,2) represent the mass (decay width) of the phys- 
ical states. The complex mixing parameters p and q obey 
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where (f> is the CP-violating phase in D°-D° mixing. 

The time-dependent amplitude for a D or D° decay- 
ing into a final state / after a time t is 



(m\ D0 (t)) 
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where A f = (f\H\D°) and A f = (f\H\D°) are the decay 
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Here A/ is the relative phase between Af and Af, and 
rd is the magnitude of the ratio between the two am- 
plitudes. Mixing will occur either if x or y is non zero, 
while CP violation in mixing is signaled by p ^ q, which 
can occur either if \q/p\ ^= 1 (CP violation in mixing) or 
4> 7^ (CP violation in the interference between mixing 
and decay). Direct CP violation is signaled by Af ^ A jr. 

In the SM, D°-D° mixing arises from |AC| = 2 (C is 
the charm quantum number) short-range box diagrams 
(see Fig. [TXcft) containing down-type quarks, strongly 
suppressed either by small 6-quark couplings (CKMsup- 
pressed) or by the GIM cancellation mechanism [T3] for 
the light d- and s-quarks. As a consequence, non-zero 
values for x and y are generated in the SM only at second 
order in SU(3)f breaking, x, y ~ sin 2 9qx {jn 2 — TO d)/ m ci 
where 9c is the Cabibbo angle and m s , m^, and m c 
are quark masses |15j . Lowest order calculations yield 
x ~ 0(1CT 5 ) and y ~ O(10~ 7 ), although enhance- 
ments due to higher orders in operator product expan- 
sion (OPE) up to 0(1(T 3 ) have been calculated [16]. 
Models involving NP can greatly increase estimates for 
both x and y [TT], but so can |AC| = 1 long-range 
SM processes with intermediate states accessible to D° 
and D° [TH] (see Fig. [TTRightV While most studies find 
| x |, | y | < 10 -3 , some estimates for x and y allow for 
values as large as O(10" 2 ) and suggest they are of op- 
posite sign [19]. Overall, theoretical predictions for x 
and y within the SM span several orders of magnitude, 
reflecting the fact that these processes are difficult to 
calculate [13II20] . 

However, it would be a sign of NP if x were to be signif- 
icantly larger than y or ifCPV either in mixing (j> q) 
or decay were observed [17j with current data samples. 
The observation of (large) CP violation as an unam- 
biguous sign for NP is due to the fact that all quarks 
building up the hadronic states in weak decays of charm 
mesons belong to the first two generations. Since the 
2x2 Cabibbo quark-mixing matrix is real, no CP vio- 
lation is possible at tree level, and only penguin or box 
diagrams induced by virtual fo-quarks can generate CP- 
violating amplitudes. However, as stated above, their 
contributions are strongly CKM suppressed. 
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Fig. 1. SM processes contributing to D°-D° mixing: (Left) 
Short-range box diagram and (Right) long-range interactions 
with intermediate states. 



3. Experimental methods 

A generic D°-D° mixing analysis is performed in three 
steps. First, the D° (or D°) flavor at production time 
(t = 0)is identified ("tagged") using D*+ -> D°tt+ de- 
cays [2T]. These events are usually selected and char- 
acterized using the invariant mass of the exclusively re- 
constructed D° meson, m D o, and the mass difference 
between the reconstructed D* + and D° mesons, Am = 
m(D* + ) — mjjo . The distribution of Am shows a narrow 
peak, due to the small Q-value of the D* + — > D°ir + de- 
cay. Other tools used to improve the event selection and 



reduce backgrounds are particle identification (for lep- 
tons, kaons or pions, depending on the D° final state) 
and cuts on D° (high) and soft pion (low) momentum. 
The charge of the soft pion from the D* decay unam- 
biguously identifies the D° flavor at production. Then, 
the D° flavor at decay time is identified using the charge 
of the final state particles. For example, if the recon- 
structed final state is a positive kaon and a negative 
pion, and the soft pion from the D* decay has a negative 
charge, then we have a "Right sign (RS)" combination, 

D° — ► K + ir~ . On the contrary, if the D meson has been 
tagged as D°, then we have a "Wrong sign (WS)" com- 
bination, D° — > K + ir~ . When we have a tagged D° me- 
son at production, a WS combination can occur either if 
the D° meson decays via a double-Cabibbo-suppressed 
(DCS) transition or if it oscillates into a D° meson fol- 
lowed by decay through a Cabibbo favored (CF) tran- 
sition. Only a time-dependent analysis of the WS rate 
allows to distinguish between these two effects. Finally, 
the production and decay vertices of the D° meson are 
reconstructed in order to calculate the decay flight length 
and hence the decay proper time t and its uncertainty 
(T t - At B factories, restricting the production point to 
the luminous region of the collider (beam spot) greatly 
improves the precision on the decay time reconstruction, 
as well as on Am. At these facilities the average decay 
length is about 240 /im, with typical resolution about 
100 ^m (the latter depends on the specific reconstructed 
final state). Analyses usually apply quality cuts on the 
proper-time error, in order to reduce effects from wrongly 
reconstructed vertices. 

The tagging can also be performed using coherent D°- 
D° production at charm factories running slightly above 
the D°-D° threshold, although in this case only time- 
integrated mixing related measurements are possible at 
present facilities [52] . On the other hand, some analyses 
can be performed using tagged or untagged samples, as 
it is the case of the lifetime differences between decays to 
CP cigenstatcs (like D° — > K + K~ ) and to the CP-mixed 
state D° — > K~n + , as discussed later. 

Four experimental techniques have been used to mea- 
sure D°-D° mixing, depending on the specific D° fi- 
nal state: WS semileptonic decays, WS hadronic de- 
cays, decays to CP eigenstates, and self-conjugate three- 
body final states containing a combination of quasi-two 
body flavor and CP eigenstates, particularly A^7r + 7r _ . 
Quantum-correlated final states at charm factories are 
also sensitive to D°-D° mixing (mostly y) via time- 
integrated observables [22] , although their sensitivity 
is not competitive with time-dependent measurements. 
These are however fundamental to provide information 
on magnitudes and phases of relevant amplitude ratios, 
as described later. 

4. Wrong sign hadronic decays 

Sensitivity to mixing using WS hadronic decays, for 
example D° — > K + ir~, is obtained by analyzing their 
proper-time evolution. Time-dependent studies allow 
separation of the direct DCS D° — > K + tt~ amplitude 
from the mixing contribution followed by the CF decay, 
D° -> D° -> K+tt- [121 [13]. Taking |A/| « 1 and 
assuming small mixing, 
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where / represents the WS final state, Rd = is the 
ratio of DCS to CF decay rates, R M = {x 2 + y 2 )/2 
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is the mixing rate, and y' = y cos Sf — xsinSf, where 
Sf = — A/ is the relative strong phase between the DCS 
and CF decay amplitudes. The minus sign originates 
from the sign of V us relative to V c d, where V denotes the 
quark- mixing CKM matrix. In Eq. ([7]), the first term 
corresponds to the DCS contribution to the WS rate 
(time independent), the second term is the contribution 
from mixing, and the third term is the interference be- 
tween mixing and CF decays. Since x and y are small, 
0(1O -2 ), it is precisely the interference term (linear in 
decay time, x and y) which gives the best sensitivity to 
mixing through y' . Let us note that Rm = (x' 2 + y' 2 )/2, 
with x 1 = x cos Sf + ysmSf. However, a direct extrac- 
tion of x and y from Eq. Q is not possible due to the 
unknown relative phase Sf. 

Searches for charm mixing in WS D° — > K + tt~ decays 
have been performed by the experiments E971 (using 
2 x 10 10 events from ir~ N interactions at 500 GeV) and 
FOCUS (from 10 6 D ->■ Kn(ir) events from 7 iV inter- 



actions), and by CLEO from 9 fb of e' 1 



T(45) 



data [23]. However, the available statistics from these 
experiments was not enough to obtain evidence of mix- 
ing. _ 

The first evidence for D°-D° mixing in WS D° — > 
K + ir~ decays has been reported by BaBar using 384 
fb -1 of data [5]. The simultaneous fit to the RS and WS 
data samples to describe the signal and the random soft 
pion and misreconstructed D° background components 
yields 1141500 ± 1200 and 4030 ± 90 signal events, re- 
spectively. Thus the fraction of WS decays is measured 
to be R ws = [0.353 ± 0.008(stat.) ± 0.004(syst.)]%. In 
the presence of mixing, Rws > Rd, as can easily be 
obtained integrating Eq. 

The measured proper-time distribution for the WS 
signal is modeled by Eq. ([7]) convolved with a reso- 
lution function determined using the RS proper-time 
fit. The proper-time distribution for WS data in the 
m D o — Am signal box is shown in Fig. [23 together 
with the fit results with and without mixing, shown as 
the overlaid curves. The mixing parameters are y' = 
[0.97 ± 0.44(stat.) ± 0.31(syst.)]% and x' 2 = [-0.022 ± 
0.030(stat.) ± 0.021(syst.)]%, and a correlation between 
them of —0.95. The ratio of DCS to CF decay rates 
is measured to be [0.303 ± 0.016(stat.) ± 0.010(syst.)]%. 
The systematic uncertainties are dominated by the sig- 
nal resolution function as extracted from the RS sample. 
As expected, Rws > Rd, revealing the presence of mix- 
ing. As another cross-check of the mixing signal, Rws 
can also be measured in slices of proper time, repeating 
the fit to the RS and WS data samples in each of these 
slices. The fitted and expected WS fractions are shown 
in Fig. [23 and are seen to increase quadratically with 
time, as expected according to Eq. (J7j). The significance 
of the mixing signal is equivalent to 3.9cr or 9.6 x 10~ 5 
confidence level (CL), where a denotes one standard de- 
viation. Separate proper-time fits to D° and D° events 
allow to determine a CP-violating asymmetry Ajj = 
(R^-Rn)/{R}^+R D ) = [-2.1±5.2(stat.)±1.5(syst.)]%, 
where is the ratio of DCS and CF decay rates 

for D a (D°), thus no evidence for CP violation is ob- 
served. 

These results have been confirmed by the CDF ex- 
periment using a data sample of 1.5 fb -1 of pp interac- 
tions at y/s = 1.96 TeV [7J. The analysis is similar to 
that from BaBar, although the different production en- 
vironment makes the details to differ significantly. The 
time-integrated fit to the RS and WS data samples yield 
(3.044±0.002) x 10 6 and (12.7±0.3) x 10 3 signal events, 
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Fig. 2. BaBar WS D° -> K+tt~ analysis 0. (a) Projections of 
the proper-time distribution of WS candidates and fit results 
allowing (solid curve) and not allowing (dashed curve) mixing, 
(b) The points represent the difference between the data and 
the no-mixing fit, and the curve shows the difference between 
the fits with and without mixing. 
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Fig. 3. BaBar WS D° — > K + w~ analysis [5]. The measured 
WS fraction in slices of measured proper time (points). The 
dashed line shows the expected WS fraction as determined 
from the mixing fit. In the absence of mixing, no time depen- 
dence would be observed. 



respectively. The ratio of WS to RS decays as a func- 
tion of the decay proper time in the range between 0.75 
and 10 D° lifetimes shows again an approximately lin- 
ear dependence, as observed in Fig. [4] The parabolic fit 
of the data in this figure returns y' = (0.85 ± 0.76)%, 
x' 2 = (-0.012 ± 0.035)%, and R D = (0.304 ± 0.055)%, 
where the errors include statistical and systematic uncer- 
tainties. The significance of the mixing signal is equiva- 
lent to 3.8er (1.5 x 10~ 4 CL). These results are essentially 
identical to those obtained by BaBar, in spite of the very 
different production environment and sources of system- 
atic uncertainties. 

An earlier search by Belle for mixing in this decay 
mode using 400 fb -1 of data did not yield clear evidence 
for mixing [24 . The time-integrated fit to the RS and 
WS data samples returns very similar yields as those 
obtained by BaBar, 1073993 ± 1108 and 4024 ± 88 sig- 
nal events, respectively, from which Rws = [0.377 ± 
0.008(stat.) ± 0.005(syst.)]%. The time-dependent mix- 
ing fit yields y' = (0.06g;g)%, x' 2 = (0.018§;gJ)%, and 
Rd = (0.364 ± 0.017)%, where the errors include statis- 
tical and systematic uncertainties. The correlation be- 
tween y' and x' 2 is —0.91. The no-mixing hypothesis 



4 



0.01r 



0.008 



0.006- 



0.004 



0.002t 



so 



and S- 



8 10 



Fig. 4. CDF WS D° -» K+n~ analysis 0. The measured ratio 
of WS to RS decays as a function of the decay proper time 
(points). The dashed curve is the result of the mixing fit, 
while the dotted line is a fit assuming no mixing. 



is excluded at 2.1a (3.9% CL). This result agrees with 
those obtained by BaBar and CDF at 2a level. Separate 
proper-time fits to D° and D° events show no evidence 
for CP violation. 

Quantum-correlated D°-D° pairs produced in 
'J (3770) decays at charm factories, with definite charge- 
conjugation eigenvalue C = — i, can be exploited to 
make a determination of the relative stro ng phase Sk-k to 
translate the measurement of y' into y |25| . At slightly 
higher energies (above DD* threshold) one can also 
produce such pairs with C = +1 (additional photons in 
the final state). One can use the fact that heavy- meson 
pairs produced in the decays of heavy-quarkonium states 
have the property that the two mesons are in CP- or 
flavor-correlated states [22]. For instance, one may tag 
one of the neutral- D mesons as a CP eigenstate through 
its decay into CP eigenstates, such Ks(n°, p°, cj, 77, 0), 
K + K~, and it + tt~. The other neutral- -D meson must 
then have opposite CP if C(D°D°) = -1 and the same 
CP if C(D°D°) = +1. Then one measures its decay rate 



into K 



which includes, as discussed, an interference 



between CF and DCS amplitudes. The measured rate 
thus depends on the CF and DCS rates and the relative 
strong phase <5jor- More generally, one can measure 
time- integrated yields of correlated ( "double tags" ) and 
uncorrclated ("single tags") neutral- D meson decays 
to CP eigenstates (CP-even and CP-odd) and flavor 
eigenstates (semileptonic and hadronic decays, with 
leptons and/or kaons as final state particles). The ratio 
of correlated and uncorrelated decay rates depends 
on mixing parameters Rd, \JRd cos Sf, y, x 2 , and 
V Rqx sin 6 f . From 818 pb _1 of data recorded at 
"J (3770) (and at slightly higher energies) and using 
branching ratios from other experiments, CLEOc ob- 
tains cosSk* = 1.03lg;?7(stat.) ± 0.06(syst.) [25]. The 
extraction of other mixing parameters is not competitive 
with time-dependent methods. 

Further evidence for D°-D° mixing has been reported 
by BaBar using a time-dependent Dalitz plot analysis 



of the multi-body WS decay D" — > A' H 



rug. The 



T(D° ^ f)(s , S+ ,t) 
e -rt 

\A f \* + \A f \ 2 (rtf^ + (Tt)\A f \\A f \y', 
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where Af(so, s+) is the DCS amplitude, Af[so, s + ) is the 
CF amplitude, and y' = y cos <5/(sn, s + )—x sin 5f(so, s+), 
with Sf(so,s + ) = &vg[A*f(so, s+)Af(so, s+)] the relative 
strong phase between the DCS and CF amplitudes, now 
varying with the Dalitz plot position. As it can be 
seen in Eq. ([5]), the sensitivity to mixing comes from 
the variation of the Dalitz plot distribution with time 
produced by the CF-mixing interference term, which in 
turn mainly depends on the interference between the CF 
D° -► K+p- and DCS D° -> K* + tt~ amplitudes, since 
these decays dominate the RS and WS Dalitz plots, re- 
spectively. 

BaBar determined the CF amplitude Af in a time- 
integrated Dalitz plot analysis of the RS decay sam- 
ple, consisting of 658,986 events with a purity of 99%. 
This amplitude is then used in the analysis of the WS 
sample, containing 3009 events with a purity of 50%, 
where the DCS amplitude Af is extracted along with 

the mixing parameters. Each of the amplitudes Af 
and Af are in turn described as a coherent sum of 
amplitudes, each describing a separate resonance (the 
usual isobar approach). Figure [5] shows the RS and 
WS proper-time distributions as well as the projections 
on so and s + . Since for both Af and Af one com- 
plex amplitude must be fixed arbitrarily and the CF 
and DCS Dalitz plots are different, the sensitivity to 
x and y is in the form y' = y cos 5 k-k-r ~ x sm &k-k-k° 
and x' = x cos Sk-ktt + V sin &Kn-K a i where 5k-k-k° is the 
strong phase difference between the DCS D° — > K + p~ 
and the CF D° — > K + p~ amplitudes. This phase is un- 
known and different from Sktt- The measured mixing 
parameters are x' = [2.61^0 68 ( stat -) ±0.39(syst.)]% and 
y' = [-0.06;^|(stat.)±0.34(syst.)]%, and a correlation 
between them of —0.75. The significance of the mixing 
signal is equivalent to 3.2cr (0.1% CL). No evidence for 
CP violation is seen. 

5. Hadronic decays to CP eigenstates 

The lifetime difference between states of different CP 
content, for example D° — ► K + K~ (CP even) compared 
to D° — > A'~7r + (CP mixed), can also be used to measure 
D°-D° mixing. For small mixing and taking A/ w 1 (for 
CP-even decays, since 6f — 0, ir for CP= +1.— 1) [12l 



(9) 



where r± = T(l + y'±) is the D° or D° effective lifetime, 

with y'± « y cos =p x sin </>. For untagged neutral-D 
mesons, 



analysis in such decays is formally similar to the WS 
D° — > K + tt~, but now the decay rate is a function of 
both the decay proper time and the Dalitz plot variables 



T (£ or D° - f C p) (t) 



(x e 



-<r±)t 



(10) 



where (r±) = T(l + (y'±)) is the average D° and D° 
effective lifetime. We clearly observe that for y ^ 0, the 
lifetimes to CP eigenstates (r±, (T±)) and CP mixed 
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purities of 98%, 99%, 92%, are shown in Fig. E](a,b,c), 
respectively, together with the projection of the fit. The 
CL of the no-mixing hypothesis (yep = 0) is 6 x 10~ 4 , 
which corresponds to a significance of 3.2a\ The mix- 
ing effect can be seen in Fig. EJd) , which shows the ra- 
tio of decay-time distributions for D° — > K + K~, ir + ir~ 
and D° — > K ~tt + decays, incre asin g linearly with time, 
as expected according to Eq. (flO]) to first order in y, 

r (D° or D° -> f CP ) (t) cx 1 - Yycpt. 




m| v [GeV 2 /c 4 ] [GeV 2 /c 4 ] 



Fig. 5. BaBar WS D° -» K+-k~-k° analysis [10]. Proper-time 
distribution for (a) RS and (b) WS decays with the fit pro- 
jection overlaid. Projections of the Dalitz plot distributions 
for (c) so and (d) s_|_. The gray histograms represent the ran- 
dom soft pion background, while the dark histograms show 
the misreconstructed D° background. 



states (r) differ. The experimentally defined observables 
are 

Vcp = 7 r-l = 7 ■ -— -lwycos^, 

\Th+h-) ( r ++ T -)/2 

A T = — — —ttxsinc/), (11) 
t + + r_ 

where t k -^+ = t = 1/T and (T h + h -) = (t + +t_)/2 with 
t± = l/r± is the mean lifetime for neutral- D mesons 
decaying into CP eigenstates. The observable A T is 
the asymmetry in their lifetimes, sometimes replaced by 
AF = tA t / (t±) = (1 — ycp)A T . In the limit of vanish- 
ing CP violation yep = y, and A T (or AY") is zero. Both 
yep an d A T (or Ay) vanish if there is no D°-D mixing. 
The measurement of yep requires precise determinations 
of lifetimes using either tagged or untagged neutral- .D 
mesons, but A T can only be measured using tagged D° 
and D° mesons. The advantage of these observables is 
that most of the systematic uncertainties related to the 
signal cancel in the ratios, although background related 
systematic uncertainties do not. 

Searches for D°-D° mixing in hadronic decays to 
CP eigenstates have been done by the E971, FOCUS, 
and CLEO experiments [25]. However, the first evi- 
dence for charm mixing in these decays has been pre- 
sented by Belle [6J simultaneously with the BaBar WS 
D° — > K + ir~ [5] evidence, and both constitute the chief 
analyses reporting the first evidences for D°-D° mixing 
(quickly confirmed by CDF [7]). 

Using 540 fb _1 of data, Belle has measured yep = 
[1.31 ± 0.32(stat.) ± 0.25(syst.)]% employing the K+K~ 
and 7r + 7r + final states (both CP-even), and found no 
evidence for CP violation in these decays since they 
obtained for the lifetime asymmetry A T = [0.01 ± 
0.30(stat.)±0.15(syst.)]% [5J. The proper-time distribu- 
tion for the D° — > K + K~, K ~tt + , 7r + 7r~ samples, con- 
sisting of 111 x 10 3 , 1.22 x 10 6 , 49 x 10 3 signal events with 
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Fig. 6. Belle D° — > K~*~ K ,n~*~n analysis [6]. Decay time dis- 
tributions of (a) D° -> K+K~, (b) D° -> K~n+, and (c) 
D° — * 7r+7r _ decays. The curves are the projections of the 
mixing fit, and the cross-hatched areas represent the back- 
ground contribution, (d) Ratio of decay-time distributions 
between D° -> K+K~, tt+tt" and D° -* R--K+ . The solid 
line is a linear fit to the data points. 



A very similar study has also been performed by 
BaBar using 384 fb _1 of data, yielding yep = 
[1.24 ± 0.39(stat.) ± 0.13(syst.)]% and Ay = [-0.26 ± 
0.36(stat.)±0.08(syst.)]% [8]. These results are obtained 
fitting simultaneously the proper-time distributions for 
the D* tagged K + K~, tt + it~, and K~tt + samples, con- 
sisting of 69, 696, 30, 679 and 730, 880 signal events with 
purities of 99.6%, 98.0%, 99.9%. The significance of 
the no-mixing hypothesis is of 3. Oct, reflecting the sig- 
nificance of the difference of lifetimes between K + K~ , 
7r + 7r~ and K~tt + , as summarized in Fig. [71 

Recently, BaBar has presented an untagged analysis 
using the same data sample of 384 fb _1 [5]. The proper- 
time distributions for the untagged D° — > K~ir + and 
D° — > K + K~ samples are shown in Fig. [H3 together 
with the projection of the simultaneous fit to these sam- 
ples. In this case, the samples contain 2710.2 x 10 3 
and 263.6 x 10 3 signal events with purities of 94.2% 
and 80.9%, for D° -> K~ir+ and D° -> K+K~ re- 
spectively. The measured lifetimes are {tk+k-) = 
405.85±1.0(stat.) fs and t k - w + = 410.39±0.38(stat.) fs, 
yielding y C p = [1.12 ± 0.26(stat.) ± 0.22(syst.)]%, which 
excludes no-mixing with 3.3o\ In this analysis, since the 
initial flavor of the decaying D° does not need to be 
identified, no D* + reconstruction is required, increasing 
significantly the reconstruction efficiency but increasing 
the amount of background. To minimize it, the life- 
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Fig. 7. BaBar D° —* K + K~ ,tt + tt~ analysis [§]. Summary of 
the measured lifetimes in the 5 D* tagged K^K~ , 7r+7r~ and 
A"~7r+ samples. 



time fit is performed in a narrow D° mass region around 
the nominal D° mass. The proper-time distribution for 
the main background component (combinatorial) is es- 
timated from sideband D mass regions, while for the 
small admixture of misreconstructed charm decays it is 
obtained from the simulation. Combining the tagged 
and untagged results taking into account both statis- 
tical and systematic uncertainties, BaBar finds yep = 
[1.16 ± 0.22(stat.) ± 0.18(syst.)]%. Summing statistical 
and systematic uncertainties in quadrature, the signifi- 
cance of this measurement is 4.1a. 
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Fig 



8. BaBar D° -» K+K~ untagged analysis [9]. (Left) 
D° -» K~tt+ and (Right) D° — > K~K+ decay-time distribu- 
tions with the data (points) and the projection of the lifetime 
fit overlaid. The gray and black distributions represent the 
combinatorial and charm background contributions. 



Yet another lifetime-difference analysis uses the 
KqK + K~ final state, where the 4>K% region and its 
sidebands are examined to extract CP-odd and CP- 
even amplitudes. This is in effect a measurement of 
the lifetime in the CP-even and CP-odd parts of the 
K S K + K~ Dalitz plot. Using 673 fb" 1 , Belle has mea- 
sured y C p = [0.11±0.61(stat.)±0.52(syst.)]% [27]. This 
is done by measuring the mean lifetime ton in the (f>Kg 
region (mainly CP-odd) and the mean lifetime tqff in 
the sidebands (mainly CP-even), along with the corre- 
sponding fractions foN and /off of CP-even events in 
these regions. The lifetime asymmetry in these regions 
can then be related to ycPi 



tqff — ton 
toff + ton 



Vcp 



foN — foFF 



1 + ycp(l - foN - /OFF) 



(12) 



from which relation the latter is then determined. The 
main systematic uncertainties in this analysis come from 



ON- OFF differences in the proper-time resolution func- 
tion and the selection criteria, while the uncertainty from 
the Dalitz model assumptions needed to evaluate the 
CP-even content is negligible (0.01%). 

Figure [53 summarizes all the available yep results. 
The combination of all these measurements is p erfo rmed 
by the Heavy Flavor Averaging Group (HFAG) [28] , and 
yields yep = (1.107±0.217)%, which differs significantly 
(about 5er) from zero. The combined lifetime asymmetry 
is A T = (—0.123 ± 0.248), thus there is no evidence for 
CPV. 
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Fig. 9. Summary of yep measurements. The average 
y CP = (1.107 ± 0.217)% differs significantly from zero gH]. 



6. D° -» AIItt+tt- decays 

The last employed technique to study D°-D° mixing 
involves the multi-body final state Kgir + ir~ . As in the 
case of the WS D° — > A + 7r~7r° analysis, the decay rate is 
a function of both the Dalitz plot variables s+ = m 2 K0 ff+ 

and s_ = mr K a and the D° decay proper time. With 
the usual approximations, 

r(D° ->/)(«+, a_,t) 



-17 



(13) 



\A f \ 2 (ny 



l R M + \A f \\y*-x 2 )/2 



(rt)\A f \\A f \y' : 



with y' = ycos[(S/(s+, s_) + </>] — xsin[(5/(s+, s_) + </>], 
where <5/(s+, s_) = arg[A^.(s+, s-)Af(s+, s_)] is the rel- 
ative strong phase between the D° and D° decay ampli- 
tudes to the same final state / = A!j7r + 7r~. Here, and 
contrary to the WS D° — * A" + 7r~7r° case, the strong 
phase 6f is fixed by the fact that the D° and D° Dalitz 
plots are identical (the s+ and s_ axes are just inter- 
changed), assuming CP is not violated in the D° decay. 
Thus the analysis is free of unknown phases, providing 
an unique method to simultaneously measure the inter- 
fering D° and D° amplitudes, the mixing parameters x 
and y (without rotations, and also their signs), and even 
the CP-violating parameters <f> and \q/p\. 
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The CLEO experiment pioneered this analysis us- 
ing only 9 fb _1 of data [55], obtaining the constraints 
(-4.5 < x < 9.3)% and (-6.4 < y < 3.6)% at 95% 
CL. Using 60 times more data, Belle has also per- 
formed this analysis |30j . first assuming CP conserva- 
tion and subsequently allowing for CP violation. The 
amplitudes Af and Af are described using a coher- 
ent sum of 18 amplitudes, dominated by CF D° — > 
K*-ir+, DCS D° -> K*+ir- and CP D° -> K a s p a de- 
cays. Assuming negligible CP violation, x = [0.80 ± 
0.29(stat.)+°;°7(syst.)+°;i°(model)]% and y = [0.33 ± 

0.24(stat.)^ ' 12 (syst.)^ [ 08 (model)]%. This corresponds 
to a significa nce of 2.2a from the no- mixing hypoth- 
esis. Figure 110.1 shows both the statistical-only and 
overall contours for both the CPV-allowed and the CP- 
conservation cases. No evidence for CP violation is 
found. 



no CPV (stat. only 
no CPV 

CPV (stat. only) 
CPV 




x (%) 



Fig. 10. Belle D° -> K%ir+ir~ analysis [30]. 95% CL contours 
for (x,y): dotted (solid) corresponds to statistical (statisti- 
cal and systematic) contour for no CPV, and dash-dotted 
(dashed) corresponds to statistical (statistical and system- 
atic) contours for the CPV-allowed case. The point is the 
best fit result for the no CPV case. 



7. Wrong sign semileptonic decays 

The most straightforward although not the most sensi- 
tive way to search for charm mixing is to use WS semilep- 
tonic decays, for instance D° -> K^+i'Pi [3T] l32l [33] . 
In this case WS combinations can only occur through 
mixing, 



T(D° -> f)(t) « e- rt (Tt) 2 R 



(14) 



Therefore, these final states are only sensitive to Rm ~ 
O(10 -4 ). Using semileptonic decays for mixing searches 
involves the measurement of the time-dependent or time- 
integrated rate for the WS decays. The main experi- 
mental challenge in these analyses is the limited mass 
resolution on Am due to the presence of neutrinos. Sig- 
nificant improvements on Am resolution are obtained 
applying kinematic constraints on the invariant mass of 
the neutrino and the kaon-lepton-neutrino system. Th e 
best current limits are from BaBar, Rm < 0.12% [32] , 
and Belle, R M < 0.06% [33], both at 90% CL, using 344 
fb" 1 and 492 fb _1 of data, respectively. 

8. Combined results 

The task of combining the wide var iety of c harm mix- 
ing results is done by the HFAG [28] . Figure 111. I shows 



the (x, y) contours of the collective experimental data, 
for the case of CP conservation. The central values 
x = (0.989 ± 0.241)% and y = (0.809 ± 0.160)% ex- 
clude the no-mixing point with 10.2o\ Other relevant 
combined parameters arc Rd = (0.3360 ± 0.0084)%, 
S Kn = 0.44 ± 0.17 rad, and 8 Kn „o = 0.24 ± 0.37 rad. 
When CP violation is allowed, the mixing parameters 
remain basically unchanged, x — (0.976 ± 0.249)%, 
y = (0.833 ± 0.160)%, R D = (0.3367 ± 0.0086)%, S K7r = 
0.46 ± 0.17 rad, and 8 K ^o = 0.26 ± 0.37 rad, and the 
following values for the CP-violating parameters are ob- 
tained: \q/p\ = 0.866 ± 0.160, <\> = -0.148 ± 0.126 rad, 
and A D = (-2.2 ±2.4)%. 
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Fig. 11. Two-dimensional contours for the mixing parameters 
(x, y) from the HFAG combination 1281 using all available 
measurements, for the case of CP conservation. The no-mix- 
ing point is excluded with a significance clearly exceeding 5cr. 



9. Summary and conclusions 

More than thirty years after the discovery of the D° 
meson [34] and the first theoretical discussion on mix- 
ing and CPV in the charm sector [4], BaBar, Belle and 
CDF Collaborations have provided compelling experi- 
mental evidence for D°-D° mixing. Collective experi- 
mental data favor the mixing hypothesis at 10.2cr level 
(including systematic uncertainties). The mixing mea- 
surement from D° lifetime differences yep is significantly 
positive (about 5cr), indicating that the \D\) eigenstate 
(ft CP-even) has a shorter lifetime than the I-D2) eigen- 
state (ft CP-odd). 

However, no observation (more than 5c) in a sin- 
gle measurement has yet been presented. In addition, 
to date there is only one direct measurement (from 
D° — > A'^7r + 7r~ decays) of x and y free of rotations 
(thus the only indication of the x and y relative sign), 
which do not differ significantly from zero but do affect 
the combination of all mixing measurements due to the 
large uncertainties that arise from the unknown phases 
Sf that are inherent in other determinations. This is es- 
pecially true for x. Clearly, more such measurements are 
a high priority and are foreseen in the future. 

The measured values of the mixing parameters x ~ 
y « 1% are about compatible with SM expectations, 
although with large theoretical uncertainties. There is 
no evidence for CP violation in D° mesons, either in 
mixing, in decay or in interference. 

Significant improvements in precision are foreseen in 
the short term with the analysis of complete data sets 
from current facilities (B factories and Tevatron). In the 
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long term, facilities about to start or just starting (LHCb 
and BESIII) have the potential to improve the precision 
on the mixing parameters in about a factor 5 and look 
deeper into CP violation searches. In the longer term, 
SuperB factories could be the last opportunity to observe 
CP violation in D°-D° mixing [35] and search for NP in 
FCNC in the down-quark sector. 
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